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INTRODUCTION 


The better we are able to predict airloads , the more efficient the 
airframes we are able to design. Our ability to predict airloads is im- 
proved through the analysis of flight measurements and wind tunnel tests , 
and through the development of theoretical analysis. 

The air vehicle development process has several phases: (1) the ana- 

lytical investigation of a large number of candidate configurations in 
the preliminary design phase; (2) wind tunnel tests to fine-tune the 
configuration and to establish its rigid aerodynamic characteristics; (3) 
analytical predictions of the effects of flexibility; (4) design, fab- 
rication, and ground test of the prototype vehicle; and (5) flight test, 
which is the final proof of the concept. 

Currently, large computer-aided design systems are being developed 
with the objective of automating as much of the development cycle as 
possible. These systems are completely dependent on the accuracy of the 
analytical methods used. Realistic evaluation of the accuracy of these 
methods is obtained from comparison with experimental results, especially 
flight test. 

Since the instriunentation, calibration, and flight testing of full 
scale air vehicles are expensive, most aeronautical research is conducted 
through wind tunnel testing and computer analysis. It is, however, 
necessary to obtain flight measurements from time to time to evaluate 
these prediction techniques. 

Presented herein are flight test results for the space shuttle 
orbiter (OV-101) wing loads measured during the approach and landing test 
(ALT) program conducted at the NASA Dryden Flight Research Center (DFRC) 
at Edwards, Calif, during 1977. The flight test results are compared 
with both wind tunnel predictions and the predictions from a FLEXSTAB 
analysis computer model. From these comparisons, an assessment is made 
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regarding the current capability for airload measurement and prediction 
on a low aspect ratio double delta aerodynamic configuration. 

Recognition is due to Rodney Rocha , of the NASA Johnson Space Center 
(JSC), for supplying flight and wind tunnel data; Donald Black, DFRC, for 
programming support; Kenneth Iliff, DFRC, for supplying basic aircraft 
parameters; and to the Air Force Flight Test Center (AFFTC) for the ground 
tracking data. 
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SYMBOLS AND ABBREVIATIONS 
load equation coefficients 


wing airload bending, N-m (in- lb) 

wing net bending, N-ra (in- lb) 

vehicle reference span, m (in.) 
wing reference span, m (in.) 


wing bending load parameter, 



lift coefficient, ^ 
qS 

pitching moment coefficient. 



A 

wing torque load parameter, 

i 

^A 

wing shear load parameter, — 

q 

vehicle reference chord, m (in.) 
wing reference chord, m (in.) 

fuselage station 
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2 2 

acceleration due to gravity, m/sec (ft/sec ) 
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lb 


HM 


ob 
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V. 


V 

W 

w 

w 

w.s. 

X, Y 


hinge moment of inboard elevon, N-m (in-lb) 
hinge moment of outboard elevon, N-m (in-lb) 
altitude, m (ft) 

components of moment of inertia about body axes 
kg-m^ (slug-ft^) 

total lift, N (lb) 

vehicle body length 

total pitching moment, N-m (in- lb) 
lateral load factor, g 

normal load factor, g 

2 

dynamic pressure, N/m (psf) 

2 

computed pitch acceleration, deg/sec 

2 2 

vehicle reference area, m (ft ) 

2 2 

wing reference area, m (ft ) 

wing airload torque, N-m (in-lb) 

wing net torque, N-m (in- lb) 

wing airload shear, N (lb) 

wing net shear, N (lb) 

true velocity, m/sec (ft/sec) 
vehicle weight, kg (lb) 
wing weight, kg (lb) 

wing station 

longitudinal and lateral axes , respectively 



X 

x/c 

y 

y/b 

z 

a 
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Subscripts : 


i 


a 
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e 


longitudinal or chordwise coordinate, m (in,) 

wing chordwise center of pressure location, percent 

lateral or spanwise coordinate, m (in.) 

wing spanwise center of pressure location, percent 

vertical coordinate, ra (in.) 

angle of attack, deg 

5 * 5 

equivalent aileron position, ■ . It . ? : * . ? . -f . lr ^ . ^ 

body flap position, deg 

elevon position, deg 

d d 

equivalent elevator position, ^ 

rudder position, deg 
speed brake position, deg 

net output from strain gage i, pm/m (pin/in.) 


gage number 

per degree angle of attack 
per degree of elevator position 


0 


Subscripts applied 

eg 

cp 

le 


zero angle of attack 
to vehicle coordinates: 
center of gravity 
center of pressure 

leading edge of wing reference chord 



ref 


wing reference point 


w wing center of gravity 

VEHICLE DESCRIPTION 


The orbiter configuration (fig- 1) is a deep boxy fuselage sup- 
ported on a double delta wing. Elevens are used for pitch and roll 
control. Split rudders are used for yaw control and also act as speed 
brakes. A body flap is used for additional pitch trim. Control surface 
maximum displacements are listed in table I. The reference geometry used 
for aerodynamic derivatives are listed in table 2. The mass properties 
for the ALT configurations are listed in table 3. 

The structure is aluminum and is protected against aerodynamic 
heating by a layer of insulating ceramic tiles. The wing structure (fig. 2) 
consists of a main box and a glove separated by a wheel well. The main 
box is composed of 6 spanwise spars with corrugated webs and 12 chordwise 
ribs with truss webs. The main covers are aluminum skin-stringer aft of 
the wheel well (F.S. 1191 to 1365) and honeycomb sandwich adjacent to the 
wheel well (F.S. 1040 to 1191). The glove is composed of truss stub-spars 
covered with aluminum skin. 


PREDICTION TECHNIQUES 

Airload predictions were generated by using both theoretical and 
wind tunnel results. The procedure in the development of the predicted 
results involved four steps , as follows : 

(1) Aerodynamic coefficients and stability derivatives were determined 
from theoretical calculations or wind tunnel force data. 

(2) Surface pressure distribution and resulting shear, bending, and 
torque loads were determined for given values of angle of attack, 
Mach number, and surface control position from theoretical 
calculations or wind tunnel pressure surveys. 

(3) Angle of attack and control surface position were determined from 

a simulator for a given flight condition using results from step (1). 

(4) Final airloads were computed for the flight condition using 
results from steps (2) and (3). 
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FLEXSTAB Analysis 

The FLEXSTAB computer program (ref. 1) provides a potential 
flow, finite element aeroelastic solution for arbitrary shapes 
at subsonic and supersonic Mach numbers. The solution produces surface 
pressures, airframe deflections, and stability derivatives. The stability 
results have been found to be quite satisfactory for conventional high 
aspect ratio configurations (ref. 2) but somewhat less satisfactory for 
low aspect ratio configurations like the YF-12 (ref. 3) and the orbiter 
(ref. 4) because of the nonpotential flow for these configurations. 

The FLEXSTAB orbiter model (fig. 3), which is similar to the model 
developed in reference 4, consists of 20 control points along the body, 

80 panels on the wing, 42 panels on the vertical tail, and 9 panels for 
the body flap. The model shown was used for subsonic Mach numbers only. 

Above Mach 1, the body flap was deleted, and for Mach 1.5 and above, the 
nose was made more conical (surface slope was decreased) to avoid mathematical 
difficulties. The model was initially run without camber or thickness for 
the lifting surfaces (that is, as a flat plate). Later camber or thickness 
was introduced by using the airfoil geometry from references 4 and 5. This 
produced marked improvement in trim, hinge moment, and wing chordwise center 
of pressure predictions. Since reference 4 found little aeroelastic effect, 
only the rigid case was examined. To predict loads, the surface pressures 
were integrated by using a follow-on computer program to obtain shear, bend- 
ing, and torque load components at the flight measurement stations. Solu- 
tions for several a , 6^ combinations were integrated to obtain wing loads . 


Then coefficients were generated to express the loads as functions of a 
and 6^. For example. 


A Bq Bg e 

e 


T = (C + C a + C„ a^)q 
^ ^0 ® 


Elevon hinge moments were also expressed in this fashion. These equa- 
tions permit the prediction of the loads for any flight condition where 
q, a, and 6^ are defined. Body flap, rudder, and speed brake positions 

were assumed to be zero for the analysis of the selected pitch maneuvers. 
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Wind Tunnel Data 


Models from 0.004 to 0.0405 scale were tested in tunnels at 
NASA Langley, Ames, and Johnson, at the Arnold Engineering Development 
Center (AEDC), and at Rockwell (Los Angeles Division) to obtain force 
and surface pressure data (refs. 6 and 7). The surface pressures were 
edited, adjusted to agree with force data, and integrated to obtain 
wing loads for various a, 6^ combinations. Load coefficients were derived 

from these results in the same form as used for the FLEXSTAB results. 


Trim Predictions 

To obtain predictions of a and 6 for the dynamic flight test cases, 

stability derivatives were entered into a standard six-degree-of-freedom 

simulation computer program which yielded time histories of a, 6 , and N 

e’ z 

for a given pilot input time history. The ALT control stick steering con- 
trol mode (CSS) was represented in the simulation. 


FLIGHT TEST TECHNIQUE 


Instrumentation 

The parameters available from the onboard recording system tape for 
defining flight condition were Mach number, angle of attack, control 
surface positions, cockpit acceleration, and rate gyro output. Ground 
tracking provided altitude, velocity, and dynamic pressure. 

Calibrated strain gages were used to measure wing loads. As shown 
in figure 4, a line of gages was installed along W.S. 134, with axial 
gages on the upper and lower spar caps , and rosettes on the spar webs and 
upper and lower skin covers. Elevon actuators were instrumented with 
pressure transducers to determine hinge moment. 


Calibration 

In accordance with the procedure presented in reference 8, which was 
used successfully on the YF-12 airplane (ref. 3), a series of concentrated 
loads was applied to the instrumented wing at the locations indicated in 
figure 5. The resulting gage outputs were recorded on a sensitive data 
acquisition system and used to derive a set of linear equations expressing 
the wing shear, bending, and torque loads (referenced to W.S. 134 and 
F.S. 1135) due to air and inertia loads acting outboard of W.S. 134. The 
derivation followed the least-square multiple linear regression technique 
described in reference 8. 


The load equations take the form: 
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Flight Data Processing 


The data from the onboard recorder were reduced to engineering units 

at JSC and sent to DFRC on tape. After correction for zero shift, the 

strains were substituted into the load equations from the calibration to 

obtain loads. Pitch rate was differentiated to obtain pitching acceleration. 

Pitching acceleration and cockpit linear acceleration were combined to 

obtain load factor (N ) at the center of gravity. Finally, the net loads 

21 

were corrected for inertia loads to obtain airloads using the following 
relations: 


V. = V.. + (N - 1)W 
A N z w 


B. = + (N. - - y„.) 


w ^ref 


^A = 


T-t + (N 
N 2 


1)W (x - X \p) 
w w ref 


where the wing weight and center of gravity are as given in table 4. 
Eleven hinge moments were computed from measured pressures and eleven 
geometry. These data, along with condition parameters, were output to 
printer and plotter. 


Flight Maneuvers 

The ALT program consisted of a series of flights where the orbiter, 
mounted on top of a modified B-747 airplane, was carried to an altitude of 
6096 to 7620 meters (20,000 to 25,000 feet) and released for unpowered 
glide flights terminating in landings on Rogers Dry Lake at Edwards , 

Calif. ALT flights were designated inert captive (IC), captive active (CA) , 
or free flight (FF), depending on the orbiter 's status (inert or active) 
and whether the orbiter was actually released (captive or free). Flights 
FFl to FF3 were flown with a tail cone attached to reduce the buffeting of 
the B-747 empennage and were thus not completely representative of the or- 
biter flight configuration. However, flights FF4 and FF5 were flown with 
tail cone off, corresponding to the configuration to be used for actual 
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orbital flights. FF5 was chosen for this study. Figure 6 shows the time 
history of altitude, speed, and dynamic pressure for flight FF5 from lauhch 
to landing, taken from ground tracking data. The launch was made at an al- 
titude of 5944 meters (19,500 feet) and Mach 0.55 and lasted for 1.8 minutes. 

Figures 7(a) and 7(b) show selected parameter time histories for the 
flight. In figure 7(a), the upper graph plots angle of attack and Mach 
number, the middle graph plots body flap and speed brake position, and 
the lower graph plots elevator and aileron position. From the middle 
plot it is apparent that the body flap moves from -11.7° to 0.6° at 
launch and that the speed brake was opened about halfway through the 
flight and closed just prior to touchdown. The effect of the speed brake 
activity can be seen in the elevator position time history and in the 
dynamic pressure plots of figure 6. In figure 7(b), the upper and middle 
graphs show rudder position and side load factor, while the lower plot 
shows normal load factor and pitching acceleration. The times for longi- 
tudinal and lateral inputs were selected from these plots. 


Load Maneuvers 

For FF5, a pitch maneuver was performed, as noted in figure 7(b), 

2 

at Mach 0.52 and a dynamic pressure of 11,060 N/m (231 psf). It reached 
a maximum load factor of 2.47g. Expanded time history plots for the FF5 
pitch maneuver are shown in figures 8(a) and 8(b). An interval was chosen 
for plotting a, 8 , and wing loads versus N from these pitch time his- 

6 Z 

tories for comparison with the wind tunnel and theoretical predictions 
(fig. 8(b)). 


RESULTS AND DISCUSSION 


Comparisons of predicted and measured results were made for aerody- 
namic derivatives, angle of attack, elevon position, elevon hinge moment, 
and wing airload and center of pressure. Predictions were obtained from 
FLEXSTAB theory (corresponding to the preliminary design mode) and from 
wind tunnel data. 


Derivatives 


Comparisons of theoretical and wind tunnel derivatives are shown 
in figure 9. The angle of attack derivatives (C- , C^ ) show the best 

a 

overall agreement, peaking at Mach 1.2. Derivatives with respect to 


elevator position (C 


L ’ 


) were overpredicted near Mach 1 but were 
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relatively close elsewhere. Aerodynamic coefficients (G^ , ) showed 

^0 “o 

the poorest correlation, with predicted C- values biased positively and C„ 

^0 "o 

predictions exhibiting a discontinuity near Mach 1. The lack of correla- 
tion in the transonic range is not surprising, considering the linear 
potential flow aerodynamics in FLEXSTAB. 


Trim 

Angle of attack and elevator position were obtained from simulator 
computer runs by using wind tunnel derivatives , theoretically predicted 
derivatives, and the recorded pilot pitch control input. Time histories 
for the FF5 pitch maneuvers are shown in figure 10. Flight measured values 
are compared with wind tixnnel and theoretical predictions. The wind tunnel 
values closely follow the flight measurements , while the theoretical values 
underpredict load factor and elevator position. 

Angle of attack and elevator position are plotted against load factor 
in figure 11. The figure shows that the slopes of a/g and S^/g are predicted 

quite well by both wind tunnel and the theoretical methods. 


Eleven Hinge Moment 

Wind tunnel and theoretically derived values are plotted with flight 
measurements in figure 12. The wind tunnel predictions compare closely 
with the flight measurements for the inboard surface and are somewhat low 
for the outboard surface. The theoretical predictions are in close 
agreement with the flight measurements for both surfaces. 


Wing Load and Center of Pressure 

Shear (V^) , bending (B^) , and torque (T^) loads were computed as 

described above. These_values were used to compute the spanwise 
(y/b) and chordwise (x/c) locations of the airload center of pressure as 
shown in figure 13. The variation of these wing load parameters with 
load factor are shown in figure 14. The three sets of data are for 
flight, wind tunnel, and theoretical results. 

The load varies linearly with load factor. The spanwise center of 
pressure location is relatively constant, while the chordwise center of 
pressure is a strong ftmction of elevator position, moving ahead of the 
leading edge (x/c < 0) for negative (trailing edge up) elevator positions. 

Wind tunnel predictions give a load/g slope close to the flight 
measured value but overpredict the value by a nearly constant increment. 


Wind tunnel predictions of spanwise center of pressure are somewhat out- 
board but agree closely on chordwise location. 

FLEXSTAB predictions are in good agreement for load/g slope. Theoreti- 
cal predictions of center of pressure location are somewhat inboard spanwise 
but close chordwise. 


CONCLUDING REMARKS 


The results of this study led to the following observations regarding 
the current capability for airload measurement and prediction on a low aspect 
ratio double delta aerodynamic configuration. First, strain gage load 
measurement was successful for a low aspect ratio wing. Second, simulation 
using both wind tunnel and theoretical derivatives produced useful predic- 
tions of load factor, angle of attack, and elevator position. And third, 
FLEXSTAB surface pressure predictions yielded useful predictions for 
both wing loads and elevon hinge moments . 


Dryden Flight Research Center 

National Aeronautics and Space Administration 
Edwards, California, April 1981 
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TABLE 1. - CONTROL SURFACE LIMITS 


Eleven deflection, positive trailing edge down, deg . . . . . . . , -35 to 20 
Rudder deflection, positive trailing edge left, deg .......... ±22.8 

Speed brake deflection, positive trailing edge open, deg ..... 0 to 87.2 


Body flap deflection, positive trailing edge down, deg . . . . -11.7 to 22.55 


TABLE 2. - REFERENCE GEOMETRY 

S, (ft^) ... . . 249.91 (2690) 

b, ra (in.) . . ... 23.79 (936.68) 

c, m (in.) ..... 12.06 (474.81) 
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TABLE 3. - ORBITER 101 FF5 MASS PROPERTIES 

W, kg (lb) 68,422.6 (150,846) 

m (in.) . ..... 27.76 (1092.8) 

z , m (in.) 9.45 (372.0) 

I^, kg-m^ (slug-ft^) . . . 1.007 X lO^ (0.743 x lO^) 

kg-m^ (slug-ft^) . . . 7.510 x lO^ (5.539 x lO^) 

kg-m^ (slug-ft^) . . . 7.827 X 10^ (5.773 x lO^) 

1^8-'"^ (slug-ft^) . . . 0.175 X 10^ (0.129 x lO^) 

TABLE 4. - ORBITER WING REFERENCE PROPERTIES 


W^, kg (lb) ....... 5825.94 (12,844) 

y^, m (in.) 6.0 (236.2) 

x^, m (in.) 31.74 (1249.8) 

c^, ra (in.) 9.21 (365) 

b^, m (in.) \ . 8.48 (334) 

Xj^g, m (in.) ...... . 28.32 (1115) 

x^ef. ® (in-) 28.83 (1135) 

^ref’ ® • 2.67 (105) 
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Front view 


Figure 3. FLEXSTAB aerodynamic model. 
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Time of day 

Figure 7. FF5 flight time history. 
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